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Abstract 



A formalism for describing an all-sky map of the polarization of the cos- 
mic microwave background is presented. The polarization pattern on the sky 
can be decomposed into two geometrically distinct components. One of these 
components is not coupled to density inhomogeneities. A non-zero amplitude 
for this component of polarization can only be caused by tensor or vector met- 
ric perturbations. This allows unambiguous identification of long- wavelength 
gravity waves or large-scale vortical flows at the time of last scattering. 
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With the COBE detection of large-angle anisotropy in the cosmic microwave background 
(CMB), results from numerous balloon-borne and ground-based experiments, and the ad- 
vent of a new generation of satellite missions, the CMB is becoming an increasingly precise 
probe of the early Universe. CMB anisotropies will help determine whether density per- 
turbations (scalar modes) are the result of inflation, topological defects, or perhaps some 
other mechanism. Detection of a stochastic gravity-wave background (tensor modes) [|TJ or 
vortical motions in the primeval fluid (vector modes) would help to discriminate between 
these models. Inflation damps out vector modes but will produce some tensor modes and 
also predicts a specific relationship between the spectrum of the scalar and tensor fluctu- 
ations 0. In contrast, topological defects will produce a mixture of scalar, vector, and 
tensor modes. Scalar modes give rise to both the observed large-scale structure and CMB 
fluctuations, while in the foreseeable future we can only expect to observe the consequences 
of tensor or vector modes through their effects on the CMB. Without a model of primordial 
fluctuations, the contribution of scalar, vector, and tensor modes to the CMB temperature 
anisotropy are indistinguishable. 

However, any mechanism which produces temperature anisotropies will invariably lead to 
non-zero polarization as well |3|-|5| . As we demonstrate in this Letter, this polarization signal 
can be used to discriminate between scalar and vector or tensor metric perturbations. COBE 
has already mapped the polarization pattern with an angular resolution of 7° (although the 
data has not been analyzed), and MAP || will measure the polarization with a resolution 
of around 0.3°. 

In prior work, the auto- and cross-correlations between the Stokes parameters Q and 
U and the temperature T have been considered. Here Q and U are defined with respect 
to particular orthogonal axes on the celestial sphere. While this formalism does provide a 
complete description of the polarization, there is no rotationally invariant way to lay down 
orthogonal basis vectors on a sphere, so the meaning of QQ, QU, UU, QT, or UT will 
depend on absolute positions of the points being correlated rather than just the relative 
position. Calculations of this sort have been done with a small-angle approximation, since 
rotational non-invariance disappears when considering only a small patch of the sky. How- 
ever, this formalism is not optimal for describing the complete temperature and polarization 
correlations present in full-sky maps. 

Here we present a rotationally covariant formalism for describing the polarization pattern 
on a full sky. The Stokes parameters, defined by the 2x2 correlation matrix of the electric 
field of incoming photons, can be described as a tensor field on the celestial sphere. The Q 
and U parameters, describing linear polarization, are just given by the symmetric trace-free 
(STF) part of this tensor. For example, in spherical polar coordinates (#,</>), where the 
spherical metric is g a b = diag(l, sin 2 9), the polarization tensor is 



The "afe" are the tensor indices, and we use standard tensor notation throughout. It is 
natural to decompose the linear-polarization pattern into STF tensor spherical harmonics 
fTjJH; which constitute a complete orthonormal set of rank-2 STF tensors on the sphere. 




(1) 



There are two types of harmonic STF tensors, Y,f m)ab and Y', 



(lm)ab i 



one of each for every 
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one of the usual spherical harmonics Y(/ m ) with I > 2. Two sets of tensor harmonics are 
required as there are two modes of linear polarization, Q and U. Since Compton scattering 
can produce no net circular polarization, the CMB is expected to have V = 0, and the V 
Stokes parameter will not be considered further. 

The harmonic expansion of an all-sky map of the CMB temperature and polarization 
can be written 



= l + ^aJ m) Y" (M (n), 
J o lm 



V nh (h 



T I,,, 



The mode amplitudes are given by 



X! [ a (L)^(>n)ab(fi) + a (L)^(7m.)a&(n) • (2) 



^ -^/^nT(n)y ( I m) (n), 



a (lm) — 7f 

af lm) = ^JdnV ab (n)Y^*(h), 
af lm) = ^JdhV ab (h)Y^*(h), (3) 

which can be derived from the orthonormality properties 

/ d n Y (lm)(n) Y (l'm>)(n) =S U >6 mm >, 

J rfnr ( f4 a6 (n)r ( ^(n) = 8 lv 8 mm ,, 
J dhY^ )ab (h)Y^)(n) = 5 lv 5 mm ,, 

J dhY^ )ab (h)Y^)(n) = 0. (4) 

Here To is the cosmological mean CMB temperature and we are assuming Q and U are 
measured in brightness temperature units rather than flux units. 
The two geometrically distinct tensor harmonics are 

Y(lm,)ab — Nl (Y(i m y Mb - -g a bY(i m y c c ^j , 

(Y(i m y ac € C b + Y(i m y. bc € C a^ . (5) 



(lm)ab c\ 



Here iVj = y2(l — 2)1/ (I + 2)! is a normalization factor, e ab is the completely antisymmetric 
tensor, and ":" indicates a covariant derivative on the sphere. In two dimensions, any STF 
tensor can be uniquely decomposed into a part of the form A, ab — {l/2)g ab A, c c and another 
part of the form B :ac e c b + B :bc e c a where A and B are two scalar functions. This decomposition 
is quite similar to the decomposition of a vector field into a part which is the gradient of a 
scalar field and a part which is the curl of a vector field; hence we use the notation G for 
"gradient" and C for "curl". Since the Y(; m )'s provide a complete basis for scalar functions on 
the sphere, the Y9, ab s and Yn m \ ab s provide a complete basis for G-type and C-type STF 
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tensors, respectively. This G/C decomposition is also known as the scalar/pseudo-scalar 
decomposition ||. 

In (9, 4>) coordinates, where Eq. (JI]) holds, the harmonics are given explicitly by 
Y G (f l \ = ^i( W ( lm ) X {lm) sine \ 

(Jm)«6 W 2 ^ X(im) sin q _ W(im) sin 2 q J , 

y c (*\ = W( ~ X ( 1 ™) W {lm) sme\ 

(Zm)«6 W 2 I W(im) gin Q X(im) sin 2 q I W 



with the definitions 



d m 
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2im ( d A 
X ^"sin^(M- COt T (;m) - (7) 

The exchange symmetry {Q, U} <-> {[/, — Q} as G^C indicates that V^)^ and V/^ a6 
represent polarizations rotated by 45°. 

A most useful property of the G/C decomposition is that, in linear theory, scalar per- 
turbations can produce only G-type polarization and not C-type polarization. This is in 
contrast to tensor or vector metric perturbations which will produce a mixture of both 
types. To understand why scalar metric perturbations do not produce a C-type polariza- 
tion pattern, consider a scalar perturbation with single Fourier mode k in the z direction. 
The polarization in a given direction can be represented by a mag nitude V = (Q 2 + U 2 ) 1/2 
and an orientation angle a from the axis defining the Stokes parameters (here, choose 9), 
where tan 2a = U/Q. For scalar perturbations, the orientation of the polarization can be 
determined only by the direction of k: thus a = if the polarization orientation is along the 
direction of k, or a = it/ 2 if the orientation is perpendicular to the direction of k. In either 
case, in a given region of the sky all of the orientations are parallel and thus the polarization 
pattern has no curl. Since the curl is a linear operator, summing over Fourier modes does 
not alter this conclusion. For tensor and vector perturbations, the azimuthal symmetry in 
the scalar case is explicitly broken, and thus the Fourier vector does not completely define 
the direction of the polarization orientation. Another way to state this argument is that 
scalar perturbations have no handedness so they cannot produce any "curl" , whereas vector 
and tensor perturbations do have a handedness and therefore can. 

Finding a non-zero component of C-type polarization in the CMB would provide com- 
pelling evidence for significant contribution of either vector or tensor perturbations at the 
time of last scattering. Given a polarization map of even a small part of the sky one could in 
principle test for vector or tensor contribution by computing the combination of derivatives 
of the polarization given by P afe : & c e c a which will be non-zero only for C-type polarization, i.e. 
when vector or tensor perturbations are present. Similarly only G-type polarization con- 
tributes tO V ab :ab - Of course, taking derivatives of noisy data is problematic; more robust 
measures are given below. 

We now turn to statistics of CMB polarization. If the cosmological inhomogeneities are 
Gaussian random noise, then to the extent linear theory is valid, the CMB fluctuations 
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will also be Gaussian random noise. Regardless of whether the distribution is Gaussian, 
rotational invariance requires that the 2-point correlations be of the form 



T \ 
a (l'm>) ] 


= Cf Su'5 mm ' , 


a (Z'm')) 


— C l 5u'S mm ', 


a fl'm')) 


— C t 5u'S mm f, 


a fl'm')) 


= Cj G 5u'S mrn 


a (£'m')) 


= Cj c 5w5 mm 


a («'m')) 


— l^l OlliO mm 



a (lm) 

G* 
a (lm) 

C* 
a (lm) 

a (lm) 

a (lm) 

a (im) a («'m')) = C GIJ 5u<5 mm i . (8) 

If we also require that the distribution of inhomogeneities be invariant under parity, then 
Cfc = Cp c = since the Y(i m ) and the Y G m ^ ab have parity (—1)' while the K^L ab have parity 
(— Measuring a non-zero Cj c and/or Cp c would be quite interesting, indicating a 
handedness to the inhomogeneities in our universe. However we do not expect this and will 
henceforth only consider the four angular power spectra {C; T , Cp, Cp, Cf G }. The first is the 
well-known angular power spectrum of temperature anisotropies while the last three, new 
to this paper, are related to various quantities in previous work (see Ref. 0). Note that 
the scalar, vector, and tensor contribution to each of the C;'s adds in quadrature, i.e. for 
X — X, G, C, TG 

£jX ^(Xscalar _|_ ^(Xvcctor _|_ ^jXtensor 

and this is true whether or not the fluctuations are Gaussian. We have argued that C ; Cscalar = 
0. 

Given an all-sky temperature/polarization map, one can determine the a^ m ^s using 
Eq. (||), and then construct estimators for the CVs in the usual way, i.e., 

l \ n T 12 l I G 12 

C T = ^ \ a (lm)\ _ \ a {lm)\ 



. Z |„C |2 i T* G 

c c _ j2 {lm) c? G = lm) } ( 10 ) 

m=-l 2/ + 1 m= _ ; 21 + 1 

If only part of the sky is mapped, the same techniques developed to analyze anisotropy with 



incomplete sky coverage |L0j may be applied to polarization to construct other estimators 



of the various CVs. The mean square polarization is 



Q2 + jj2 = 2V ab V ab = Vl + V 2 c (11) 

where 



V 2 °° 91 4- 1 - — - V 2 °° 9/ 4- 1 - 

J 1=2 S7F J 1=2 87F 



Since scalar modes do not contribute to Vq, this statistic provides a powerful and unam- 
biguous model-independent probe of tensor and vector perturbations. 
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To test a given spectrum of tensor modes against a polarization map, comparing the 
complete set of predicted Cp with the estimators Cp is more powerful than considering only 
Vq, if the detection has sufficient signal-to- noise. Usually, however, the theory being tested 
has scalar as well as non-scalar modes, along with undetermined cosmological parameters. 
If so, the most information can be extracted from the map by comparing the entire set of 
predicted moments, {Cj , Cp, Cp, Cp G }, with the measured estimators |TT|P^| . 

Note that only the C\ . s potentially allow detection of a small vector or tensor signal. If 
scalar perturbations dominate, then the vector or tensor signal in {Cp, Cp, Cp G } may be 
swamped by the cosmic variance in the scalar signal, but the Cp's are not contaminated in 
this way. The cross-correlation moments Cp , which differ for scalar, vector, and tensor per- 



turbations [13], will be larger than the polarization autocorrelation moments. Furthermore, 
the noise in the temperature and polarization maps will be to a large extent uncorrelated, 
so the pixel-noise variance to these moments will be substantially reduced. Therefore, the 
temperature-polarization cross correlation may be measured with some precision. 

Much of the small-angle formalism of Refs. [p|-p|,p!3]-[T6| can be reproduced by replacing 
the Y(z m )(n)'s in our formalism with Fourier modes, e lln , and using regular derivatives rather 
than covariant ones. This small-angle formalism is completely analogous to that developed 
above and will provide an accurate description of a region of sky small enough to be approx- 
imated by a flat surface. The G/C decomposition in the small-angle formalism can be used 
to detect non-scalar perturbations on small scales, though the tensor and vector signal are 
liable to drop off rapidly at angular scales smaller than a few degrees. 

To make contact with previous work, we can write the two-point temperature and po- 
larization correlation functions [|3|-§].|l3|-[l6j in terms of multipole moments ||. Although 



correlation functions of Stokes parameters which appear in the previous literature depend 
on the positions of the points being correlated, rotationally-invariant correlation functions 
exist which are closely related to those discussed above. To construct them, define Stokes 
parameters Q r and U r with respect to axes which are parallel and perpendicular to the great 
arc (or geodesic) which connects the two points being correlated. The two-point correlation 
functions are 

C T (9) = (T(n 1 )T(n 2 )) fii . fi2=cose 
= ^oE^^(cos£), 

C«(0) = (a(fii)Qr(fi 2 )} firfi2=cose 

2/ -I- 1 

= T^^iCpW^O) +i CPX {l2) (6,0)}, 
C u (9) = (U r (ni)Ur(ri 2 )) &1 . &2=cose 

91 4- 1 

= T^^ r [cpw m (e,o) + i cpx m (e,o)i 

C^(6) = (T(nOQ r (n 2 )) fil . fl2=cose 

= T^ 2 -^N lC rP?(cos6), (13) 

where 9 is the angle separating the two points and is the m = 2 associated Legendre 
function. Since Q r and T are invariant under reflection along the great arc connecting the 



6 



two points while U r changes sign, (Q r U r ) = (U r T) = if statistical invariance under parity 
holds. Eqs. (|13|) reduce to the correct small-angle formulae [16] when 6^1. The functions 
{C T (6),C Q (6),C U (6),C TQ (6)} are a different way of representing {Cp, C ; T , Cp, C^ G } and 
vice versa. In Gaussian models either set provides a complete statistical description of the 
temperature and polarization patterns. 

The largest hurdles to detecting and characterizing CMB polarization are sensitivity and 
foregrounds. In adiabatic models with standard recombination the polarization is only a 
few percent of the anisotropy, although it may be larger in reionized ||15|| , isocurvature, or 



topological-defect models. Thus the polarization signal is at least an order of magnitude 
below current experimental sensitivities. Experiments planned or envisioned over the coming 
decade, however, will likely attain the raw sensitivity necessary for detailed polarization 
investigations. Polarized emission from foreground sources is a relatively unknown factor 
at this time. Foreground emission and any Faraday rotation will certainly contribute to 
the C-type polarization, but these contaminants can be subtracted using multi-frequency 
observations. On sub-degree scales, where the signal from vector and tensor modes are liable 
to be negligible, any measurable C-type polarization is a likely indicator of contamination. 
Polarization measurements will be difficult, but the promise of using them to detect gravity 
waves or vorticity, and hence to discriminate between cosmological models, makes these 
measurement potentially very valuable. 
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